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In this paper we study the critical behaviour of the fully-connected p-colours Potts model at the 
dynamical transition. In the framework of Mode Coupling Theory (MCT) , the time autocorrelation 
function displays a two step relaxation, with two exponents governing the approach to the plateau 
and the exit from it. Exploiting a relation between statics and equilibrium dynamics which has been 
recently introduced, we are able to compute the critical slowing down exponents at the dynamical 
transition with arbitrary precision and for any value of the number of colours p. When available, we 
compare our exact results with numerical simulations. In addition, we present a detailed study of 
the dynamical transition in the large p limit, showing that the system is not equivalent to a random 
energy model. 



I. INTRODUCTION 

Mean-field spin-glass models can be divided into two 
main classes, the ones which undergo a continuous tran- 
sition and the ones which, instead, display a jump in 
the order parameter. In systems belonging to the former 
class, at a certain temperature T s a second order phase 
transition takes place, with a continuous growth of the 
Edwards-Anderson order parameter qsA = jf J2i (Si) 2 
and zero magnetization (in absence of magnetic field): 
the spins are essentially frozen in a random direction so 
that the global mean magnetization vanishes while the 
mean squared magnetization is finite. In the low tem- 
perature phase the replica symmetry is broken with a 
continuous pattern (Full RSB) or with a step-like pat- 
tern (1-RSB) according to the Parisi scheme [1] and the 
order parameter is, in fact, a non-trivial function q(x). 
One can also study the Langevin dynamics of these sys- 
tems, showing that exactly at the thermodynamic tran- 
sition temperature T s there is ergodicity breaking, there- 
fore we can say that, in correspondence of the static tran- 
sition, a dynamical transition takes place too. 
There exists another class of mean-field spin-glass mod- 
els (like the p-spin or the p-colours Potts model [2, 3]) 
which display two different transitions: at a temperature 
T s there is a thermodynamic phase transition which is of 
the second order in terms of potentials but can be discon- 
tinuous in the EA order parameter. The low temperature 
phase is (at least in the vicinity of the critical tempera- 
ture) 1-step replica symmetry broken. At a temperature 
Td > T s a dynamical phase transition occurs, where the 
system's relaxation time becomes infinite and the ergod- 
icity is broken [2, 4]. This is due to the fact that at the 
dynamical transition the equilibrium state splits into a 
large (exponential in the system size) number of excited 
states, represented by free energy local minima. Since in 
mean-field the barriers between these states become in- 
finitely high in the thermodynamic limit, the equilibrium 
dynamics remains stuck forever in one of them and the 
overlap cannot relax to zero. 

This second class of mean-field systems has been shown 



to share some relevant properties of structural glasses [5— 
8] , more specifically, the dynamical equations are exactly 
equivalent to those predicted by the Mode Coupling The- 
ory (MCT) above the mode coupling temperature T mc 
where ergodicity breaking occurs. The analogy between 
structural glass models (with self-induced frustration) 
and proper mean field spin-glasses (with quenched dis- 
order) has been widely studied and has provided rather 
accurate predictions [9-12] . In systems with continuous 
transition, above T s the spin-spin time correlation func- 
tion C(t) = (<Tj(0)c7j(f)) decays exponentially at large 
times, which means that the system is ergodic. Lowering 
the temperature the relaxation time grows until it di- 
verges exactly at T s (the static transition temperature) 
so that the ergodicity is broken and the relaxation (at 
large times) follows a power law C(t) <~ t~ v with some 
exponent v. 

The systems belonging to the discontinuous class intro- 
duced above, behave quite differently: above Td the time 
correlation function displays at first a fast decay to a 
plateau and then a slow decay to zero (in absence of a 
magnetic field) [4]; the length of the plateau grows low- 
ering the temperature until it diverges at Td- 
According to MCT the approach to the plateau and the 
decay from it are both characterized by a power-law be- 
haviour, respectively 

C{t) ~q d + ct- a (1) 

C(t) ^q d - c't b (2) 

where qd is the height of the plateau and the two expo- 
nents satisfy the exact MCT relation 

r 2 (i-«) = r 2 (i + b) = 

r(i-2o) r(i + 26) 1 ' 

and A is usually treated as a tunable parameter (see for 
example [13]). 

The exponents a and b have been computed exactly only 
for the spherical p-spin model [4] because the dynamical 
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equations are particularly simple and correspond to the 
so called schematic MCT models. 

In most of the cases it is instead very difficult or im- 
possible to compute the exponents in a purely dynami- 
cal framework, both analytically or through Monte Carlo 
simulations. 

Numerical simulations are often difficult to interpret and 
give quite poor indication of the value of the exponents 
due to strong finite size effects; if the system is not in- 
finite, barriers between metastable states cannot be in- 
finitely high, the dynamics does not remain stuck in a sin- 
gle state and all the observables eventually relax to their 
equilibrium value since, through activated processes, the 
configuration is able to explore the whole phase space. 
The extent of this effect depends on the specific model 
we consider and, in particular, on how fast the barriers 
between metastable states grow with the size of the sys- 
tem. 

Recently, a connection between the mode coupling ex- 
ponents and some purely thermodynamic quantities has 
been introduced [14]; this connection suggests a quite 
simple recipe to compute the dynamical exponents ex- 
actly starting from the static mean-field theory which is 
much easier to work out for all the reasonable models one 
can think of. 

The aim of this paper is to apply this technique to the 
disordered Potts model and compute the MCT exponents 
for any value of the number of colours p. The Potts model 
is particularly interesting because, as will be pointed out 
in the following sections, the parameter p allows to switch 
from a continuous transition (p < 4) to a discontinuous 
transition (p > 4) , moreover in the latter case it works as 
a tuning parameter for the magnitude and separation of 
the static and dynamical transitions. In the following we 
will not make use of the symplectic representation which 
is widely exploited in literature [15-17]. The outline of 
the paper is the following: in section II we give a sketch 
of the technique used to compute the MCT exponents 
in a generic model, in section III we summarize some of 
the necessary known results about the Potts model, in 
section IV we compute the dynamical exponents for the 
Potts model for arbitrary value of the parameter p, in 
section V we compare our theoretical exact results with 
numerical simulations and, finally, in section VI we give 
our conclusions and final remarks. 



II. HOW TO COMPUTE THE EXPONENT 

Given a fully-connected model it is possible to compute 
the Gibbs free energy as a function of the order parame- 
ter which, in the case of a spin-glass transition is the well 
known overlap matrix Q. The thermodynamic value of 
the order parameter can be determined minimizing the 
Gibbs free energy functional. It can then be expanded 
around the replica symmetric saddle point solution, giv- 
ing raise to eight different kinds of third order terms. For 



our purposes, only two of them will be relevant, namely: 
Wi Tt(SQ 3 ) = Wl ^ 6Q ab 6Q bc 5Q ca (4) 

a,b,c 

and 

( 5 ) 

a.b 

In the case of continuous transitions it has been found 
[14] that there exists a quite simple relation between the 
exponent v and the two coefficients w\ and w 2 - 

r 2 (i-^) = w 2 (t s ) 

T(l-2v) Wl (T s ) [) 

In the case of discontinuous transitions it can be shown 
[14] that a relation analogous to (6) holds at the dynami- 
cal transition which, again, gives the connection between 
the dynamical exponents a and b and the static coeffi- 
cients, namely 

r 2 (i-«) = r 2 (i + b) = w 2 (T d ) 

r(l-2o) r(l + 26) Wl (T d ) ( ' 

where, differently from the former case, the expansion of 
the Gibbs free energy has to be performed around the 
dynamical overlap (the height of the infinite plateau at 
the dynamical transition). 

In order to compute the two coefficients W\ and w 2 one 
must determine the expression of the Gibbs free energy 
as a function of the overlap and then expand it to third 
order around the RS thermodynamic value q. The reason 
why the expansion has to be performed around a replica 
symmetric solution will be clarified in section IV. In fully 
connected models, introducing a replicated external field 
e, the free energy reads 

/(£) = - ^ In J dQ exp N (S[Q] + Tr eQ) (8) 

which, for N — > oo, can be evaluated at the saddle point 

f(e) = - ^cxtr Q (S[Q] + Tr eQ) (9) 

We can immediately notice that the equation above ex- 
actly defines f(e) as the Anti Legendre Transform of the 
effective action 

/( £ ) = ALT.S[Q] (10) 

and, again, by definition the Gibbs free energy T(Q) is 
the Legendre Transform of /(e), yelding 

r(Q) = LT/(e) = LT (ALT S[Q]) = S[Q] (11) 

This implies that the functional form of the Gibbs free 
energy is exactly the same of the effective action. In 
fully connected models, we can then directly expand the 
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latter. 

The general form of the third order term in the effective 
action reads 

s (3) = £ w ab 

,cd,ef SQabSQcdSQef (12) 

(ab)(cd)(ef) 

Since a ^ b, c ^ d and e 7^ / and the coefficients W 
are computed in RS ansatz, we can have eight different 
vertices: 



with 



Wap^j^a = Wl , Waji^.afj = W 2 
W a /3 tal 3 ta ~ ( = W3 , W a p ta p n s = W4 



(13) 



Following Ref. [18], Eq. (12) can be rephrased in the 
following way 



a@~f aft 
+ W 3 ^2 SQap&QaP&Qa-t + W 4 £ 5Q a/3 SQ a^SQ^g 

+ Wb $Qa05Qaj8Q{35 + Wq ^ SQ a pSQ cn 5Q a s 

(14) 



with 










tUl = 


W 1 - 


3W 5 + 3W 7 - 


W 8 




w 2 = 


\w 2 


- 3W 3 + ^W 4 


+ 3W 5 + 2W 6 


- 6W 7 + 2W S 


w 3 = 


3W 3 - 


- 3Wi - 6W 5 - 


- 3W 6 + lbW 7 


-6W 8 


= 


\ Wi 


- \w 7 + *w s 






w 5 = 


3W 5 - 


- 6W 7 + 3W 8 






w 6 = 


w e - 


3W 7 + 2W 8 






W7 = 


\w 7 








w 8 = 


lw 8 









(15) 

It is therefore sufficient to compute the eight W coeffi- 
cients and use Eq. (15) to get w\ and w 2 . 



III. THE POTTS MODEL: SUMMARY OF 
KNOWN RESULTS 

We consider the p-colours disordered Potts Hamilto- 
nian 



U = - JijVivuaj) 

<i,j> 



(16) 



77(0,6) = pS a<b - 1 



(17) 



where p is the number of colours and a = 0, 1, • • • ,p — 1. 
The sum is extended over all the possible couples taken 
from N spins and the couplings are independent gaus- 
sian random variables with mean Jq/N and variance 
J 2 /N, where the normalization is needed in order to ob- 
tain a finite thermodynamic limit. As usual, we are inter- 
ested in computing the mean-field free-energy exploiting 
the well known replica trick in order to average over the 
disorder 



InZ 



lim — In Z n 

n->0 n 



(18) 



Carrying on the computation we obtain the replicated 
partition function in a functional integral form 



''•' = 1 



DQDm cxp(-NS[m, Q}) 



(19) 



where the "effective action" S[m, Q] is a function of two 
order parameters: the magnetization to" and the overlap 
Q"s, with greek replica indices a, f3 = 1, • • • ,n and latin 
color indices r, s = 1, • • • ,p. 

5[TO, Q ] = ^(l- P ) + ^££(Q r t) 2 

a<f3 r.s 



2p 



Jo + P-T 



4 

2P 



H[m,Q,{<j}} 



+ 



(3 2 J 2 



££(TO«) 2 -lnTr w e«^M] 

a r 

(20) 

a<f} r,s 

Jo + 13 J 



2P-2 



££m«r?(a«,r) 

a r 

(21) 



In order to determine the order parameters we can use 
the two saddle point equations, which read 



m? = (( V (a a ,r))) 



(22) 
(23) 



where ((•••)) is the average taken with respect to the 
measure 

e «[m,Q,{ CT }] 

MW) = Tr {T} ettKQ.M] (24) 

The order parameters are clearly redundant, in fact they 
satisfy the following constraints: 



r 

E< = ° 



(25) 
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In the particular case p — 2 one recovers the SK model 
solution [1]. 

For p > 2 ferromagnetic ordering is always preferred be- 
low some temperature Tp [3]. An upper bound Tp for 
the temperature Tp below which ferromagnetic ordering 
appears is [19] (from now on we consider J = 1): 



2(1 - Jo) 



(26) 



For p > 4, in order to prevent that ferromagnetic order- 
ing occurs at a higher temperature than the spin-glass 
one, the couplings should be antifcrromagnetic in aver- 
age, with J less than some (negative) threshold value. A 
lower bound for the critical mean- value is Jp = (4 — p) /2 
[19]. Under this condition the magnetization is zero and 
it is straightforward to show that, as a consequence, the 
overlap has the symmetry Q"f = Q a/3 r](r,s). Then it is 
possible to write the Gibbs free energy as a function of a 
unique overlap matrix in the following way [19]: 



r[Q] = J(p-i)/5 2 E^- 



— log Tr exp 



a<fS 
a</3 



(27) 



Differentiating with respect to Q a p one obtains the sad- 
dle point equation 



Qap — 



1 Tr V (a a , a" ') exp (U[Q, a}) 
p-1 Trexp(U[Q,a}) 



(28) 



with 



U[Q,a}=p 2 J2Q a ^(a a ,a^) (29) 

a</3 

It has been shown [3] that for 2.8 < p < 4 the system un- 
dergoes a continuous transition at a temperature T s = 1 
with 1-step RSB. The breaking point is m = (p — 2)/2. 
For p > 4 the transition occours at a temperature T s > 1, 
it is discontinuous and the RSB scheme is 1-step with 
breaking parameter m — 1 at criticality. In this case 
there exists a (dynamical) glass transition, associated to 
the static one, occurring at some temperature Td greater 
than T s . The static and dynamical transition temper- 
ature and overlap can be determined numerically with 
great accuracy using the marginality condition and the 
techniques described in Rcf. [19]. We briefly summarize 
the results here. 

We can compute the free energy (27) in the 1-RSB ansatz 
with qi — q and qo = and expand it at first order around 
m = 1 as To + (m — l)T\(q), 



T(q) = -/3 2 (1 - p) - log(p) + (m - l)x 
(\/3 2 (p - l)q 2 + \l3 2 q(p + 1) + Iog(p) - h 



(30) 



where the integral I 2 is given by, 



, (3 2 pq, 
exp( ^— ) 



n 

r=l 



dy r 



2tt 



e PVm>yi log 



^2 ex p(f3(qp) l2 yr) 



vr=l 



(31) 



For m = 1 the expression (30) gives the high-temperature 
free energy which is independent of q. This general ex- 
pansion allows to determine the static and the dynamic 
transition. 

The static temperature is determined imposing that a 
solution q s exists, satisfying the following conditions: 



<9r N 



0Ti 
dq 

(ri) 



= 







(32) 
(33) 



On the other hand, for the dynamical transition temper- 
ature, we must search for a marginal stability and the 
condition becomes. 



\ dq 



o 



d 2 T 



d 2 v 1 

dq 2 



(34) 
(35) 



In the language of the Franz-Parisi potential [20] the 
two conditions above correspond, respectively, to the ap- 
pearence of a local minimum (horizontal flex) for the dy- 
namical transition, and to the fact that this minimum 
reaches the same height of the paramagnetic one, for the 
static transition. 

As we will sec in the following, the p-dimcnsional integral 
I2 in (30) is extremely hard to evaluate numerically as 
soon as p > 2. Therefore in Ref. [19] the authors use the 
identity 



log(l + A) 
and taking 

A 



dx 



'(1-e 



-Ax\ 



= expC%p) = y r ) - 1 



they obtain the result 

4 - r ± 

Jo x 

with 

dy 



x {l-e x w(xe^ qp )w p - 1 (x)} 



(36) 



(37) 



(38) 



r 

w(x) = I 

J —I 



exp 



1 



-y -xexp 



which is much easier to evaluate numerically [27] . 



(39) 
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FIG. 1: Black solid line: y = q. Dashed lines: y = L (oo) (£,<?) 
for £ = 4 , 4/3. Blue solid line: y = L (oo) (£, g) for £ = 2. 
The dynamical transition is located in £ = 2, for which the 
line y = q is tangent to the curve y = l/ 00 '^, g). 



Infinite number of colours 

It has been pointed out [3] that the Potts model be- 
comes a Random Energy Model (REM) in the limit 
p — > oo, with a critical temperature that diverges like 



T„ = 



1 



log(p) 



(40) 



In the following, we show that the limit model is not 
exactly a REM. The first of Eq.s (32), which is satisfied 
both at the dynamical and statical transition, can be 
written in the following way 



1 



p-1 



p[j ' G P {z)^2 ex P (^(P'?) 1 



/ E^ 1 exp(2/3( W ) 1 / 2 z r ) 
ELiexp(/3( M )V 2 z r ) 




(41) 



g=^L.(pL(p)( / 8,g)-l) 
P- 1 



(42) 



We have been able to show that, as a function of the 
rescaled temperature 



2 P 



log(p) 



(43) 



the right hand side of equation (41) tends to a Heaviside 
function in the infinite p limit (see Appendix), 



(44) 



with breaking point q = l/£. 
As can be easily seen from Fig. 1 both qd and q s go to 



1 in the limit p —¥ oo and the dynamical transition is lo- 
cated at the rescaled temperature such that the breaking 
point of L^ 00 ) is 1. 



Td 



21og(p) 



(45) 



While in a random energy model (REM) the ratio be- 
tween Td and T s is formally infinite [21], in the large -p 
Potts model this ratio tends to a finite value, namely 



1.414- 



(46) 



therefore, the limit model is still a "glassy" model with 
a dynamic and a static transition. 

This is at variance with the Ising p-spin model in the 
p — > oo limit that goes to a REM [21]. 



IV. THE POTTS MODEL: MCT EXPONENTS 

The determination of the mode coupling exponents fol- 
lows essentially the steps described in section II. Ex- 
panding the effective action (27) to third order around 
the replica symmetric saddle point we obtain the eight 
coefficients 



Wi 


= Ri 


-3(p 


- 1)?M 2 - 


f 2(p 


- i)V 






w 2 


= R-2 


-3(p 


- l) g M x - 


f 2(p 


- 1)V 






w 3 


= R a 


-(P- 


l)qM x - 


2(p- 


l)gM 2 - 


H2(p- 


i)V 




— i?4 


-(P- 


l)qMt - 


2(p- 


l)qM 3 - 


^2(p- 


i)V 


w 5 


= R 5 


-2(p 


- l)qM 2 - 


-(P- 


1)«M 3 - 


H2(p- 


i)V 




= Rs 


-3(p 


- l)qM 2 - 


f 2(p 


- 1)V 






w 7 


= R7 


-(P- 


1)«M 2 - 


2(P- 


l)qM 3 - 


^2(p- 


i)V 


W 8 


= i?8 


-3(p 


- 1)?M 3 - 


f 2(p 


- 1) V 







(47) 

where the replica symmetric overlap is determined 
through the saddle point equation 

«={(l(^/)» (48) 
the "mass matrix" can assume three different values 



M 1 = {{ V (a a ,^)r,{a a ,aP))) 
M 3 = (( V (a a ,^)r,(a->,a s ))) 



(49) 
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and the six replica cumulants are given by 



Hi 


1 i m ( „a 

= \\V{v 




(T>)r)(a> 


))) 


R 2 


= \\V{v 


—B\ rr ,l —a 




a ))) 


^3 


//„/„.<* 
= {{Viv 


a' J )r)[p 


a H )i]((T 


u ))) 


i?4 


= ((V(° a 


(j fi )r](a a 




<^)» 


R 5 


= M° a 


a fi )n{a fi 


ct 7 )?7(ct 7 






= (v7K 




o~<)r)(o a 


* 4 )» 


R7 


= (foK 


a%{a a 






R& 


= ((v(* a 






0» 



Given the relationship (15) between W\, w 2 and the W 
coefficients one obtains: 

W 1 =Ri- 3i? 5 + 3i? 7 - #8 

W2 = \[ R -2- 6^3 + 3i? 4 + 6R5 (51) 

+ 4(i? 6 - 3R 7 + R 8 )} 

where only the disconnected cumulants are left. 

If the thermodynamic phase transition is continuous, 

then it coincides with the dynamical one (as in the SK 



Lh = -Lkio 
V 



Lk — -i^L k oo 



withC7 p (z) = (n^ 1 ^exp(-i, r 2 )). 

As already pointed out, the above result holds both for 
continuous and discontinuous transitions with the only 
difference that in the former case q and the L integrals 



model) . In this case dynamical quantities at infinite time 
relax to their static value [22] and the averages above can 
be computed in a replica symmetric ansatz taking finally 
the limit n — > 0. If, instead, the transition is discontin- 
uous then the coefficients have to be computed at the 
dynamical transition, where quantities at infinite time 
do not relax to their equilibrium (thermodynamic) value 
but remain stuck at their value inside the most excited 
metastable states. The averages should then be com- 
puted inside a single state; this corresponds to taking 
a 1-RSB ansatz with breaking parameter m — > 1 or, if 
the mutual overlap between different states is as in our 
case, a RS ansatz with the number of replicas n — > 1 
[23] . Finally, we can assume replica symmetry and leave 
n unspecified, obtaining the expression for the two coef- 
ficients. 

w 1 =p 3 {L 3 -3U + 3L 23 -L 222 ) (52) 

and 

p 2 

w 2 = — (1 - q) + 

+ \p 3 (l - 6A3 + IOL4 + 3^22 - 12i 23 + 4L 222 ) 

where, exploiting the fact that (n(a a , a b )) 2 = (p — 1) + 
(p — 2)n(<j a ,a b ) the saddle point equation becomes 



(53) 



(54) 



I 

are computed at n — while in the latter we consider 
n = 1. 



1 



' [ g P fe) (EjU exp (P{pg)^z r )) n - 2 (gjU exp (2P(pq)^z r )) ' 
/SpUKE^iexp^) 1 /^))" 



and we have defined the class of integrals 

1 



^klh 



fG P (2-) (E^exp^V^)) 1 



g p (z) (^exp {p{pqY' 2 z r ) 



n—k—l — h 



^Texpfaipq^zM x (^cxp(?/3( M ) 1 



I 2 - 



vr=l 



vr=l 



Z r X 



J^cxp (hfiipq) 1 



/ 2 ; 



\r=l 
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The continuous transition 



The discocontinuous transition 



If p < 4 the phase transition is second-order [3] with 
q{x) continuous for p = 2 and step-like for p = 3. In the 
case of continuous transitions we have to consider n = 
and q = and the result (already found in Ref. [3]) is 
very simple, namely: 



w 2 

Wi 



P 



(55) 



which yclds v 2 = 0.5 and v z ~ 0.395. 
As in the case of the fully-connected model, it can be 
proven [24] that, on the Bethe lattice, for p < 4 the phase 
transition is second-order. The difference on the Bethe 
lattice is that, for p = 3 and low enough connectivity, 
the order parameter q{x) is a continuous (Parisi type) 
function while for high connectivity it becomes a step- 
like function (as in the fully-connected case). This does 
not affect the result which is again (for p < 4) given by 
Eq. (55). 



We are interested here in the case p > 4, when the 
system undergoes a dynamical transition, therefore we 
must take the limit n — > 1 in order to compute the 
exponents correctly. 

The computation of the overlap q and of the third order 
coefficients w\ and u>2 involve p-dimensional integrals 
(54), which become very difficult to evaluate numerically 
as soon as p is grater than 2. In order to overcome this 
issue, using the identity 



1 

Jk 



(k 



-l)Uo 



x k - x e~ Ax 



(56) 



the integrals (54) can be rewritten in following form, 
which is more suitable for numerical evaluation 



Lkih 



e 2 



5/5 -pq 



- [p e ^+i+hfp^ Pq I dx x k+i+h-2 wP -i (x) w ^ e (fc+l+fc)/J a M ) + 

vr- 2 (x)w(xe( k +W 2 ™) w (xe h ? 2 ™) + 



:x k+l + h-2 w ,. 



. e (k+h)f3 2 pq} w ( xe lp A pq 



p{k + l + h-2)V 

+ p(p-l)^ (fe+ ° 202p<? e^ 2 ^ / dx: 

Jo 

/>oo 

+ p(p-l)e5( fe + /l ) 2 ^ e ^ 2 ^ J dxx k+l+h - 2 W P- 2 (x)w(xt 

+ p(p - 1) e W+*f0 2 P1 e hk 2 2 P<l jT dx x k+l+h-2 wP -2 {x) w ^ e {l+h)P m ^ w e kf Pq ^ + 

+ p(p- l)(p- 2)e^ k2+l2+h2 ^ 2 P" J™ dxx k+l+h - 2 w p - 3 (x)w(xe k P 2p o') w (xe 1 ^) w (xe 1 ^') 



(57) 



with w{x) given in equation (39). 

Through identity (56) we have been able to reduce p- 
dimensional to "sort of" 2-dimensional integrals. 
They are not technically 2-dimensional integrals because 
in formula (57) for each value of the integration variable 
x we have to perform 2, 3 or 4 integrations to obtain the 
function w in different points (instead of just one inte- 
gration which would be needed in a regular 2-d integral) . 
Some care is needed in the computation of w(x), espe- 
cially for small x's since the integrand has an extremely 
steep growth near zero and the integration step must be 
taken very small. 

In order to go to very large values of p it should be bet- 
ter to recast the integrals (57) in their asymptotic form 
using Eq.s (A4) and (A7) given in the Appendix. We re- 
port below the results for different values of the number 
of colours p. 



p 


T d 


Qd 


A 


a 


5 


1.0101 


0.09507 


0.8764 


0.2290 


7 


1.0577 


0.2206 


0.8236 


0.2651 


10 


1.1420 


0.3238 


0.8052 


0.2759 


12 


1.1970 


0.3665 


0.8002 


0.2787 


15 


1.2748 


0.4114 


0.7962 


0.2810 


20 


1.3926 


0.4598 


0.7930 


0.2827 


30 


1.5941 


0.5142 


0.7904 


0.2841 


40 


1.7648 


0.5455 


0.7895 


0.2846 


100 


2.4964 


0.6187 


0.7892 


0.2848 



V. COMPARISON WITH NUMERICAL 
SIMULATIONS 

For p = 10 there are Monte Carlo simulations per- 
formed in Ref. [25, 26] to investigate the finite-size ef- 



8 



fects on the glass transition. They find that the ther- 
modynamic static quantities such as the energy, the en- 
tropy, the susceptibility and the overlap distribution dis- 
play very strong finite size effects. It is found also that 
the system remains always ergodic and the plateau in 
the equilibrium spin-spin correlation function C(t) is al- 
most invisible even at temperarures close to the dynami- 
cal transition Tq and big system sizes. Since for TV — > oo 
the physics of the system should be descibed by the exact 
mode-coupling equations, they expect a divergence of the 
relaxation time r(t) with a power law behaviour at the 
dynamical transition. 



T (X 



(58) 



with an exponent 7 which, in mode coupling theory is 
related to the exponents a and b through the exact rela- 
tion 



1 1 

2a~ + 2b 



(59) 



They plot t~ ~ for a set of reasonable trial values of 7 and 
find that the data are linearized in the region 1.1 <T< 
1.4 for 7 = 2.0 ± 0.5. This value of 7 gives, through the 
relation (59), an indirect estimate for a(j) w 0.36. Two 
different kinds of finite-size scaling are considered in order 
to perform extrapolations of C(t, N) at N — > 00. They 
find that only one of the two gives a C(t) which is well 
compatible with a power law behaviour of the type (1). 
In this way they can make a rough direct estimate of the 
exponent, obtaining a = 0.33 ± 0.04 which, despite the 
difficulties (identification of the plateau, extrapolation 
etc..) is close (within 2a) to our exact computation. 



VI. CONCLUSIONS AND REMARKS 

In the first part of this paper we have presented a re- 
view of some known results about the disordered Potts 
model. In the second part, exploiting a technique that 
has been recently developed [14], we have computed the 
dynamical exponents of the autocorrelation decay both 
in the case of continuous and discontinuous transition, in 
a completely static framework. 

In Fig. 2 we show the plot of the exponent a as a function 
of the parameter p from p = 5 up to p = 40. 

Since our computation is non-perturbative, the expo- 
nents can be determined with arbitrary precision and 
they can be taken as a reference in numerical simulations. 
Knowing a priori the exponents in the thermodynamic 
limit, one has an additional tool for studying, for exam- 
ple, the finite size effects and the deviations from MCT in 
a numerical simulation of a finite (fully-connected) sys- 
tem. 



0.28 



0.26 



0.24 



0.22 1 — ' — 



1 20 30 40 50 60 70 80 90 100 



FIG. 2: The exponent a computed for different values of the 
number of colours p in the discontinuous regime. 
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Appendix A: Infinite p limit 

In this Appendix we compute the dynamical critical 
temperature T<j and overlap qd in the limit p — > 00 and 
the first correction. 

Using identity (56) the p-dimensional integral L of Eq. 
(42) can be rewritten in the following way 

f°° r / \ 1 p~ 1 / \ 

L (p) (^q)= dx w(xe-^ 2pq ) w(xe^ 2pq j 

° (A1) 
Setting a = fip^^q 1 ' 2 and making the change of variables 



ka 2 +ax 



X = e 2 



(A2) 



-a -\-ax 



p-i 



(A3) 



the integral becomes 

f dxae^+^wU^^A \w(t 

J — oo 

Through standard manipulations it can be shown that 

aeW+^w (e a2+a ^ = -^e'^ F(x) (A4) 
where Fix) tends to unity in the a — > 00 limit 

/°° r "i 2 - 

dy e y - eV e'^-^ (A5) 
- 00 

Moreover we have clearly 

(e-" 2+ai ) -> 1 (A6) 



w < 



9 



The behaviour of k; p_1 and, consequently, of the integral 
L is now determined by the leading order of the first 
correction A = w — 1, in fact 



u< 



-\-OLX 



p-i 



exp [p log (1 + A)] (A7) 



We have to compute the leading behaviour of 



A 



dye 



-2V 



exp 



making the change of variables 

y = z + a — x 

we obtain the form 



(A8) 



(A9) 




1 1 



/ ' 1 
/ 1 1 



FIG. 3: The right hand side of equation (A17) for p = 

io 10 ,io lo2 ,io lo3 ,io 101 



dz e~ i^ z ~^ 2 e~ az [exp (— e az ) — 1] Solid line: £ = 2. Dashed line: £ = 2/3 

-co 



(A10) 

which in the limit of large a goes to the following form 



/oo 
dze-^-^Oi-z) 
-oo 



l e -|^+aS Erfc / ' 



(All) 



Given this correction, Eq. (A7) becomes 



v— 1 
w F ~ exp 



2 V72 



(A12) 



In order to obtain a finite result for finite x we must have: 

„2 



a 



log(p) - y + to = 



(A13) 



for some value of t wich now becomes our variable. 
Under this condition ui p_1 behaves like 6(t — x) and the 
equation for the overlap now reads 



(A14) 



Using a rescaled inverse temperature f3 2 = £ log(p)/p we 
have 



a 2 = £9k>g(p) 
and Eq. (A13) becomes 

e?--A=(fr)*-2 = o 

v 1o s(p) 

Substituting Eq. (A16) into Eq. (A14) we get 



(A15) 



(A16) 



(A17) 



Eq. (A17) can be used to obtain approximate solutions 



in the large p limit. 

From Eq. (A17) it is clear at this stage that L^°°\ given 
a rescaled inverse temperature £, is a step function in 
q with breaking point 2/£ (see Fig. 3). It then follows 
straightforwardly that in the p — » oo limit we have 



(oo) 

Id 

Joo) 



(A18) 



We are now interested in computing the first correction 
to the critical point for non-infinite p. Considering that 
1 << t « log(p), from Eq.s (A14) and (A16) , we have 
at leading order 



9 = 1 



^-2 



-t 2 /2 



-2V2 



(A19) 



+ 



we can then define 

e = 1 — q 

satisfying the two coupled equations 
2e = 2V2- 



(A20) 



e = 



t 



(A21) 



We can obtain t 2 from the second of Eq.s (A21) and plug 
it into the logarithm of the first one getting, at leading 
order 

log( £ ) = -^log(p)( A1 -2e) 2 (A22) 

At criticality the derivative of Eq. (A22) must hold as 
well, giving the second contraint necessary to determine 



10 



(p - 2e) log(p) 



(A23) 



both e and /x 

1 

e ~~ 4 

Substituting this last equation into (A22) one obtains 

1 



e 2 log(e) 



which at leading order gives 



log(p) 



(A24) 



Given e and \x we can write the dynamical overlap and 
critical temperature with the leading correction for large 
p: 



qd = 1 - 



Tl = 



\og{p) log (log(p)) 



2Iog(p) 



1 - \/2 



log (log(p)) 
log(p) 



(A27) 



e = 



log(p) log (log(p)) 



(A25) 



Substituting into (A23) and taking the leading order we 
get the other correction 



= 2y/2 



log (log(p)) 
log(p) 



(A26) 
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